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The nature of the ghost cavity in the NGC 741 group 



o 
o 

(N 
O 

Q 



Oh 
6 



> 

O 
in 



(N 

o 



Nazirah N. Jetha 1 ' 2 *, Martin J. Hardcastle 3 , Arif Babul 4 , Ewan O' Sullivan 5 , 
Trevor J. Ponman 2 , Somak Raychaudhury 2 , Jan Vrtilek 5 

1 Laboratoire AIM, CEA/DSM - CNRS - Universite Paris Diderot, DAPNIA/Service d'Astrophysique, Bat. 709, CEA-Saclay, F-91191 Gif-sur-Yvette Cedex, France 
2 School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT 

3 School of Physics, Astronomy and Mathematics, University of Hertfordshire, College Lane, Hatfield, Hertfordshire ALIO 9AB 
4 Department of Physics and Astronomy, University of Victoria, Victoria, BCVSP 5C2, Canada 
5 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 



ABSTRACT 

We discuss the effects of energy injection into the intra-group medium of the group of galaxies 
associated with NGC 741. The X-ray emission reveals a large bubble, which in the absence of 
a currently bright central radio source, may have been inflated by a previous cycle of nuclear 
activity . If the bubble is filled with a light, relativistic fluid we calculate that in expanding, it 
has done more than sufficient work to counteract the energy lost from the intra-group medium 
via radiative cooling; the bubble can provide this energy as it expands and rises. Using upper 
limits on the flux density of the plasma filling the bubble at 330 MHz and 1.4 GHz, we derive 
constraints on its electron energy distribution and magnetic field strength. We show that the 
data require the high-energy cut-off of the electron spectrum to be very low compared to the 
cut-offs seen in more typical radio sources if the fluid filling the bubble is a conventional 
relativistic plasma. This suggests that the fluid filling the bubble may not have evolved by 
expansion or synchrotron losses consistent with a dead radio source, leaving a puzzle as to 
what the origin of the bubble may be. 
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1 INTRODUCTION 

, It is well known that there is a lack of very cool gas in the cen- 
tres of galaxy group s and clusters (e.g. IPeterson et al.l I2OO1I and 
ISakelliou et alj 20021) . even in systems where the cooling time is 
significantly shorter than the Hubble time. Why this should be so 
remains a mystery, and suggests that there must be some mecha- 
nism to prevent catastrophic cooling. 

If a bubble inflated by an active galactic nucleus (AGN) is 
over-pressured with respect to the inter-galactic medium (IGM), 
then it will expand, doing work on the IGM (see for exam- 



ple iNusser et alj|2006h . Depending on how significantly the bub- 
ble is overpressured, the expansion may be trans or s uper-sonic, 
drivin g shocks into the IGM (e.g. in Centaurus A iKraft et ail 
120031) . or the expansion may be more subtl e, doing gentle P dV 
work on the gas as the bubble expands (e.g. Begelman 2001 and 
iRuszkowski & Begelman! [20021) . Once the bubble has expanded 
and reached pressure equilibrium with the external medium, it 
will begin to rise b uoyantly. As it rises, it may distort (e.g. 
IChurazov etai]|200lh . The bubble will also do further work on the 
IGM whilst it rises. Such bubbles could remain intact until long 
after the AGN has switched off, with the spectrum of the bubble 
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plasma steepening so that the emission at higher radio frequencies 
drops to below detectable rates. 

Observationally there is much evidence for this picture. A 
number of clusters show radio signatures of jets terminating in 
lobes that, in turn, are coincident, in projection with X-ray sur- 
face brightness d epressions that, on o ccasion are bounded by weak 
shocks (c.f. (e.g. Birzan et al. 2004) and references therein). The 
simplest explanation is that the inflating radio bubbles are displac- 
ing the X-ray emitting plasma. 

There are, however, cavities with no detectable radio emis- 
sion (see for exam ple Abell 1795, lEttori et alj|2002l ; Abell 2597, 
IClarkeetalJl2005l ; HCG 62. iBirzan et alj|2004) . These so-called 
'ghost cavities' are typically assumed to be features created in the 
more distant past and whose radio emission has faded over time. 
This interpretation finds some support in the fact that in any given 
cluster, the ghost cavities are typically found further away from the 
cluster centre than the cavities with detectable radio emission (c.f. 
IBirzan et ai]|2004l) . 

Whether or not a cavity is associated with radio emission, of 
course depends critically on the nature of the fluid that fills the bub- 
bles. If the bubbles are inflated by AGN, then it would be likely that 
the fluid is a relativistic plasma. However, the lack of associated ra- 
dio emission in ghost cavity systems implies that the plasma must 
be in such a state that any emission is below detectable levels. The 
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lack of observed radio emission places strong constraints on the 
nature of the radiating particles and the field. 

In this paper we investigate a centrally located ghost cavity 
in the galaxy gro up associated with NGC 741. This is an optically 
identified group (Zabludoff & Mulchaev|[l998h . consisting of ap- 
proximately 48 galaxies within 1 Mpc, with a mean redshift of 
0.019 and a velocity dispersion a — 458 ± 66 km s . The BGG, 
identified as NGC 741, and its companion galaxy NGC 742 are sep- 
arated by 45 arcsec (17 kpc), with NGC 742 lying almost exactly 
to the E of NGC 741, together with a much smaller third com- 
panion approximatel y 9 arcsec (3.4 kpc) to the NE of NGC 741 
(iBonfanti et aT]|l995h . lMulchaev & ZabludofJ ( fl998h first reported 
a significant extended X-ray halo in this group, extending to 
19.2 arcmin from the centre of NGC 741 (approximately 440 kpc in 
our adopted cosmology), with a luminosity of 5.3 x 10 41 erg s -1 
(again in our adopted cosmology). 

The organization of the rest of the paper is as follows. We 
present the X-ray data and our analysis in Section [2] In Section [3] 
we present the radio data and also comment on the dynamical state 
of the system. In Section[4]we discuss the energetics involved and 
in Section [5] we place limits on the nature of the fluid filling the 
cavity. Our conclusions are given in Section[6] Throughout this pa- 
per, we use Ho = 70 km s~ Mpc -1 , with a ACDM cosmology. 
NGC 741 is at a redshift of 0.019, and 1 arcsec=0.38 kpc. All errors 
are quoted at the l-a level. 



2 X-RAY DATA 

The Chandra data of NGC 741 (ObsID 2223) were obtained from 
the archive, and processed as part of the study of ljetha et al.l d2007l) ; 
the data were reprocessed using CIAO 3.4 and CALDB 3.3.0 fol- 
lowing the CIAO threads online to apply new calibration products 
and a new event file was created. The data were cleaned using a 3u 
clipping algorithm to remove times of high background, leaving a 
30.3 ks exposure. 

In addition, we obtained and reduced the XMM-Newton data 
(sequenc e number 01 53030701) for the source (originally pub- 
lished in lMahdavi et alj|2005h . using SAS version 7.0.0. The orig- 
inal exposure times were 8.7 ks (MOS) and 7.0 ks (PN ). The data 
were fi ltered for flares following the method described in lPratt et al.l 
(2007). and vignetting correctio ns were applied to th e data follow- 
ing the procedure described in lArnaud et alj d200ll) . After clean- 
ing, we were left with good time intervals of 5.8 ks (MOS 1), 5.9 ks 
(MOS2) and 3.8 ks (PN). 



2.1 Spatial analysis 

We smoothed the exposure corrected and binned Chandra data in 
the 0.5-5.0 keV energy range using the task CSMOOTH with a min- 
imum smoothing signal to noise ratio of 3 a (see Fig.Q]left panel). 
This process revealed an elliptical cavity to the west of the host 
galaxy (see also Vrtilek et al (in prep)). Increasing the minimum 
smoothing significance to 4 a resulted in a somewhat smaller cav- 
ity. We use an elliptical region to approximate the boundaries of the 
cavity seen with 3 a smoothing; this region has semi-major axis, 
at ~19 kpc and semi-minor axis, bt ~14.5 kpc, (marked by the 
elliptical region in Fig[TJ. We assume a prolate ellipsoid so that the 
third axis Cf, = fej. 

To check that this cavity is real and not just an artefact of adap- 
tive smoothing we firstly created an exposure corrected image in 
the 0.5-5.0 keV energy range as above and smoothed it with a fixed 



Gaussian kernel of FWHM 5 arcsec. We then used the Sherpa fit- 
ting program to fit a single 2-D /3-model to the Chandra data to 
model the large-scale emission. This single-component fit is in- 
tended to model the large-scale group emission, allowing it to be 
subtracted from the image and revealing structures in the core. As 
a single-component model can only approximate the true surface 
brightness distribution, we do not use this model to derive physi- 
cal parameters. The model centre is slightly off-centre due to the 
structure in the centre of the group. We then created an image of 
the residuals of the fit, which is shown in Fig[T] (right panel). The 
cavity can be clearly seen, and we overlay the same ellipse as in 
the left panel to guide the eye. We further created a surface bright- 
ness profile in the 0.5-5.0 keV energy range which was centred on 
NGC 741, and covered, in azimuthal bins, the region occupied by 
the cavity, together with a comparison surface brightness profile in 
the same energy range for the entire system excluding the region 
we used to extract a profile of the cavity. The opening angles of the 
wedges used for the two profiles were 80 and 280 degrees respec- 
tively, and a local background was used for background subtrac- 
tion. The resulting profiles are shown in Fig.[2](left panel) overlaid 
on each other. There is a clear 3 a surface brightness drop which 
begins at approximately 10 kpc from NGC 741 and continues to 
approximately 40 kpc, implying that that the cavity seen in the im- 
ages is real. At the centre of the cavity, the ratio, Rsb of the surface 
brightness of the cavity, to that of the undisturbed gas, is approxi- 
mately (0.4 ±0.1). Compared to other systems, where the surface 
brightness contrast s are of the order o f 20 -30% of the surrou nding 
surface brightness teirzan et alj|2004l and lRaffertv et al.ll2006l ). the 
bubbles here are of a similar surface brightness contrast. 

For completeness, we generated a surface brightness profile in 
the 0.3-2.0 keV energy range from the XMM-Newton data (shown 
in Fig [2] right panel). We use this profile to provide a model for 
the surface brightness for use in Section [231 We used the XMM- 
Newton rather than the Chandra data for this since the XMM- 
Newton data extends out to larger radii, providing a better model 
for the large scale emission. The profile was fitted by a -model 
with f3 = 0.428 ± 0.005 and core radius, (r c = 2.1 ± 0.2) kpc. 
We also smoothed the XMM-Newton data with a Gaussian beam of 
FWHM 6 arcsec (shown in Fig.[3j to see if the cavity was detectable 
in the data. Whilst the cavity is not visible in the XMM-Newton data 
due to PSF blurring and the shorter exposure time from the deficit 
seen in the Chandra data, we would expect the surface brightness 
as viewed by XMM-Newton to drop from 0.35 counts arcsec -2 to 
0.14 counts arcsec -2 . In the XMM-Newton data, we observe a sur- 
face brightness of 0.2 ± 0.5 counts arcsec -2 in the region defined 
by the cavity, consistent at the l-a level with the cavity seen in the 
Chandra data. 



2.2 Spectral analysis 

From the XMM-Newton data, we extracted a global spectrum from 
a source region of 8.8 arcmin centred on the BGG. The source re- 
gion was chosen as we detect emission at 3-a significance out to 
this radius. The X-ray backgrou nd is subtracted using the back- 
ground files of iRead & Ponmanl J2003h an d we use the double 
subtraction method of Arnaud et al, J2002I) to take into account 
both the X-ray and the particle backgrounds correctly. The spec- 
trum was fitted with a simple absorbed MeKaL model, the de- 
tails of which are given in Table Q] We find a temperature of 
1.7 ± 0.1 keV, with metal abundances fixed to 0.3 Zq. The fit has 



147 for 115 degrees of freedom. We find a luminosity for 



the group of 6.6lo.g x 10 4 



erg s which agrees well with the 
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Figure 1. Smoothed, exposure corrected Chandra images of NGC 741 in the energy range 0.5-5.0 keV as described in Section l2T1 The left panel shows the 
CSMOOTH image overlaid with 1.4 GHz VLA contours, whilst the right panel shows the Gaussian smoothed residuals to the 2-D fit. In the left panel, the cavity 
described in Section |2~T| is clearly visible, and is marked by the black ellipse. The cavity is also clearly seen in the residual image. 
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Figure 2. Chandra (left) and XMM-Newton (right) surface brightness profiles for the NGC 74 1 group. For the Chandra data, profiles were taken in wedges in 
the direction of the X-ray cavity (open points) and the undisturbed gas (solid points). There is at least a 3 <r drop in surface brightness in the region between 
the dashed lines. The XMM-Newton profile is azimufhally averaged and uses data from the Ml, M2, and PN instruments of XMM-Newton . The profile was 
fitted with a beta-model and this fit was used in Section l2~3l to model the large scale emission of the group. 



iMulchaev & Zabludoffj i 19981) ROSAT result when differing cos- 
mologies are tak en into account. Our v alue for temperature agrees 
well with that of Ma hdavi et al.l d2005l) . who use the same XMM- 
Newton data. 

We next generated azimuthally averaged temperature profiles 
from both the Chandra and XMM-Newton data. The profiles were 
again centred on the BGG, and an ellipse defining the region con- 
taining the bubble was excluded from the Chandra data as shown 
in fig.Q] F° r the Chandra data, annular regions were defined out to 
70 kpc (corresponding to 3.4 arcmin), and a blank sky background 



of the same size and placement on the chip as the source region was 
used for background subtraction. The background files were pro- 
cessed in the same way as the data (following the instructions in the 
CIAO threads) and normalized such that both background files and 
the data had the same flux in the 8-10 keV energy band. Similarly, 
for the XMM-Newton data, annuli were defined out to 8.8 arcmin 
and background regions as descr ibed earlier were used. We then 
followed the method described in Jetha et al.l J200II) for the Chan- 
dra and XMM-Newton data separately, fitting an absorbed MeKaL 
model to each annulus to deproject the spectra, to produce the de- 
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Figure 3. XMM-Newton mosaiced image, using the Ml, M2 and PN cam- 
eras of the group associated with NGC 741. The data have been smoothed 
with a Gaussian beam of FWHM 6 arcsec. 



Table 1. Details of the MeKaL model fitted to the XMM-Newton data for 
all three instruments (Ml, M2 and PN). The Ml and M2 temperatures were 
tied together, and metalicities were fixed to 0.3 Zq. The hydrogen column 
density was initially fixed to the Galactic value, but then allowed to vary. 
This did not produce a significant difference in the fits. We fitted an area out 
to 8.8 arcmin, as described in the text. 



Component 


Parameter 






Absorption 
Absorption 
Ml 


JV H (xl0 22 cm- 
iV H (xl0 22 cm~ 


2 )(Galactic) 
~ 2 ) (Fitted) 


0.044 

0.042 ±0.02 


MEKAL 


kT (keV) 
Norm 




1.7 ±0.1 

(1.69±°;f) x lO" 3 


M2 








MEKAL 


Norm 




(1.691-™?) x lO" 3 


PN 








MEKAL 


kT (keV) 
Norm 




1.6 ±0.1 

(1.06±°;°|) x lO" 3 




L x (10 41 ergs- 
X 2 (d.o.f.) 




147(115) 



projected temperature profiles shown by the open points in Fig.|4^ 
for Chandra . We do not show the XMM-Newton points here, but 
the Chandra and XMM-Newton data are in good agreement in the 
regions of overlap. The normalisations of the fitted MeKaL mod- 
els were used to obtain the electron density profile (n e (r)), since: 



4tt [Da(1 + z)Y N mek x 1.17 
1 x 10- 14 V 



(1) 



where Da is the angular size distance to the source, z is the redshift 
of the source, N me k is the normalization of the MeKaL compo- 
nent fitted to the annulus, and V is the volume of the annulus from 
where the spectrum was extracted. 

In Fig. [4] we overplot two different deprojections. The open 
points show the deprojection when we disregard the presence of 



the cavity and treat the IGM as azimuthally symmetric. The filled 
points show the deprojection when we take the presence of the cav- 
ity into account and exclude a region corresponding to the cavity 
from deprojection. As the volume excluded due to the cavity can 
be a significant fraction of the volume of the shell in this case, this 
must be taken into account when doing the deprojection. We dealt 
with the excluded volume in calculating V in Eqn.Q]by scaling the 
volume of an individual shell according to the fractional area of the 
annulus that was used to extract the spectrum. Thus, for an annulus 
towards the centre of the bubble, the volume excluded was close 
to the elliptical volume of the bubble that contributed to the shell, 
whilst for an annulus which only contained the edge of the bubble, 
the excluded volume may have overestimated the contribution of 
the bubble to the shell. 

The temperature and density profiles are shown in the top two 
panels of Fig|4] We further derived the pressure profile of the data 
(shown in the bottom left panel of Fig.|4j using 



P(r) = 2.96 x 10~ 9 n(r)T(r), 



(2) 



where the numerical factor is used to convert the temperature and 
density to cgs units, n (r) is the electron density profile, and T (r) 
is the temperature profile. To ensure that our findings are robust, 
we redid the deprojection with different regions for the spectral ex- 
traction. We found that our profiles were generally consistent with 
each other, regardless of the extraction regions chosen for the de- 
projection. 

The entropy profile shown in the bottom right panel of Fig [4] 
shows marginal signs of disturbances caused by the bubbles in the 
gas surrounding them. This suggests that the bubbles may be doing 
work on the IGM. Comparing our pro files to those of Mahdav i et al, 
we find that our profiles and those of Mah davi et alj |2005) agree 
well within the core of the gr oup and at large radi i. However, in 
the region between 20-60 kpc, Mahda vTet aljlioo^) find a slightly 
lower entropy (but consistent within the errors) than what is found 
here. This discrepancy arises because the structure in this region 
is complex and may be significant enough to affect the results of 
a ny deprojection ; thus the different method of deprojection used 
in Mah davi et all J2005) will produce slightly different results. Our 
method of radial deprojection may well average over real asym- 
metric structure, but comparison of profiles obtained via different 
methods suggests that this does not change the final values too dras- 
tically. 



2.3 Actual location of the cavity 



In Section |2TT1 we calculate that the ratio of the surface brightness at 
the centre of the cavity compared to the undisturbed surface bright- 
ness, which we termed Rsb in Section [2~71 is 0.4 ± 0.1. We ca n 
use this ratio, using a method similar to that of IWise et al.1 d2007h . 
to estimate the three-dimensional location of the cavity within the 
X-ray emitting plasma. We assume that the undisturbed IGM of 
the NGC 741 can be described by the azimuthally averaged fit to 
the XMM-Newton data, using the /3-model as fitted in Section [2~71 
The radial distance of the bubble from the centre of this spheri- 
cally symmetric density profile is 25 kpc as measured from the X- 
ray maps (FigQJ. The unknown is the distance of the cavity centre 
along the line of sight, which we denote z (taking z = to be the 
mid-plane of the spherically symmetric group). For a given z the 
measured density profile allows us to calculate the surface bright- 
ness drop produced by a cavity of the observed dimensions, at the 
observed radial distance. 

From the data, R S b = 0.4 ± 0.1 (Section|2jT]l; if the cavity 
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Figure 4. The temperature, density, pressure and entropy profiles of the hot gas obtained from the Chandra data. The open points show the data for the 
entire group, whilst the solid points show the profiles for the group with a region containing the bubble removed. The vertical dotted lines show the extent 
of the cavity as described in Section |2TI determined by the smoothed image, assuming that the centre of the cavity is located at 30 kpc from NGC 741 (see 
Sectionl23l. 



were at z — 0, then we calculate that Rsb = 0.3, which is within 
the uncertainty on our measured value. Thus it is possible that the 
cavity is in the z = plane. We find that a cavity placed 14 kpc 
away from the z = plane reproduces a ratio of 0.4. Given the 
errors on the decrease in surface brightness, we can calculate an 
error in the location of the cavity; the bubble could be in the z = 
plane, or as far away as 22 kpc from the 2 = plane. Combining 
this with the measured radial distance of 25 kpc we find that the 
cavity is located at a physical distance of (29±4) kpc from the cen- 
tre of the NGC 741 group. Taking this uncertainty into account, the 
ambient temperature and pressure of the cavity are in the range 1- 
2 keV and 3-10.4 x 10~ 12 dyne cm" 2 respectively (taken from the 
deprojection in which the cavity was excluded: Fig[4}. This range 
of temperatures and pressures has implications for both the energy 
imparted to the IGM and the time-scale over which the energy can 
be imparted, which we discuss in Section[4] 



3 RADIO SOURCES IN THE SYSTEM 

NGC 741/2 has been known to be associated wit h a relatively bright 
radio sourc e (4C 05.10) since ea r ly wo rk by Ijenkins & Scheuea 
( 19805 and iBirkinshaw & Daviej Jl985h . Radio images showed 
two bright peaks close to t he centres of the two galaxies surrounded 
by more diffuse emission. Birkinshaw & Daviesl (l985) argued that 
all the radio emission, including the bright peak close to the nucleus 
of NGC 742, was associated with NGC 741: on this interpretation 
NGC 741 would have been something akin to a classical double 
source and the compact source near NGC 742 would have been its 
hotspot. It is probably because of this interpretation that NGC 741 
is generally thought of as being the radio galaxy of the system. 
Venka tesan et alJ d 19941) used VLA observations to show that the 
component close to NGC 742 was truly compact and well aligned 
with the centre of the galaxy, and proposed instead that the diffuse 
radio emission was a head-tail or narrow-angle tail (NAT) source 
associated with NGC 742. 

To allow us to make the most definite statements possible 
about the nature of the radio emission, we have acquired multi- 
frequency VLA observations (listed in Table[2} of the system from 
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the archive and reduced them in the standard manner in AIPS. For 
display purposes we combine these with a Hubble Space Telescope 
images taken with the WFPC2 instrument on 1997 Jun 20, obtained 
from the public archive. In the images obtained from the data (e.g. 
Fig. [5j, the two compact sources seen in earlier observations are 
well aligned with the nuclei of the two galaxies. However, twin 
jets are clearly visible emerging from the nucleus of NGC 742 and 
bending westwards to merge into the larger-scale extended emis- 
sion, while there is no small-scale extension of any kind associated 
with the nucleus of NGC 741. (A third bright radio source in the 
field, a small triple whose central component has a position of RA 
01 h 56 m 28.5 s Dec. +05° 38' 24", is coincident with a very faint 
point source in the HST image and is most likely a background 
quasar.) The lack of any extensio n in the NGC 74 1 nuclear source 
is consistent with the results of IXu et al.1 {2000), who show that 
the source is unresolved even on parsec scales. In the VLA data 
the compact components have flat radio spectra between 1.5 and 
4.9 GHz, indicating that there is self-absorbed parsec-scale struc- 
ture in both. 



Table 2. Radio observations used in our analysis 



We therefore agree with Venkates an et al.l J 19941) that it is en- 
tirely plausible that the extended radio emission in the NGC 741/2 
system comes only from a narrow-angle tail source associated with 
NGC 742. T he 400 km s~ difference in the radial velocities of the 
two galaxies dMahdavi & Gellerl2004l) is consistent with a narrow- 
angle tail interpretation of the radio emission, since the much more 
massive NGC 741 is likely to be close to the centre of mass of the 
group (its radial velocity is the same as the mean of the 48 galaxies 
within 1 Mpc - 5570 km s -1 ). The common stellar envelope of the 
two systems (clearly visible in Fig. [5} implies that they are close 
enough to be interacting tidally, but they need not be as close as 
they appear in projection (18 kpc): if NGC 742 were displaced by 
a few tens of kpc along the line of sight with respect to NGC 741 
it would help to explain the lack of correspondence between the 
X-ray emission (dominated by gas on smaller scales) and the ob- 
servable radio emission, and in particular the absence of any X-ray 
cavities associated with the radio source. NGC 741 clearly hosts 
some form of active galaxy, but has no evidence for current large- 
scale jets. 

We further argue that the bubble visible in the system is as- 
sociated with a previous outburst of NGC 741 rather than with 
NGC 742 for several reasons. Firstly, NGC 742 already shows ev- 
idence of old emission in the form of a long tail; if the cavity had 
been caused by NGC 742, then the bubble must be older than the 
current large scale emission. However, the projected distance from 
the core of NGC 742 to the end of the tail is 160 kpc, while the 
bubble is located only 45 kpc away. If the bubble is from a previ- 
ous outburst of NGC 742, then we would have expected it to have 
travelled further than the observed tail from the current outburst. 
Secondly, to have got to its current position relative to NGC 742, 
the bubble would have had to have crossed the centre of the group 
against the pressure gradient. Thirdly, if the radio emission from 
NGC 742 were removed from the picture, NGC 741 would look 
much more similar to some other well- studied 'ghost cavity' sys- 
tems like NGC 4636 dJones et alj|2003) or HCG 62 dMorita et all 
2006), with a weak central radio source bearing no obvious relation 
to the observable cavity. Therefore, in the remainder of the paper 
we adopt the hypothesis that the bubble was caused by a previous 
outflow of NGC 741 and that the location of the large scale emis- 
sion of NGC 742 is entirely coincidental. 



VLA obs. 
ID 


Date 


Configuration 


Frequency 
(GHz) 


Time on source 
(min) 


AB920 


1999 Jul 18 


A 


1.5 


50 


AB593 


1991 Dec 16 


B 


1.5 


10 






C 


4.9 


26 


AH276 


1988 May 19 


C 


1.5 


30 






c 


4.9 


76 


AS827 


2005 May 07 


B 


373 





4 ENERGY IMPARTED TO THE IGM 

In Section Q] we discussed the possibility that bubbles will do 
mechanical work on the IGM by expanding and rising, pre- 
venting catastrophic cooling. Here, we calculate the amount of 
energy available to do this. The physics of A GN-inflated bub- 
bles is discussed in length in t he l iterature (e.g. Begelmar]|200l[ 
iRuszkowski & Begelmanll2002l and lNusser et alj|2006l) . and in all 
cases, the initial overpressure factor is unknown. We now briefly 
investigate the total amount of energy availabl e based on whethe r 
the bubble was inflated gently (as discussed in lAllenetal .1 (2006), 
where the bubble is overpressured by a factor of 2 with respect to 
the IGM), or if the bubble was injected more violently, say if the 
overpres sure factor was 10 , as seen in some young low-power radio 
sources dKraft et al.ll2003h . 

Assuming that the bubble lies 29 kpc from the centre of 
NGC 741 and that the ambient pressure is in the range 3 — 10 x 
10" 



dyne cm 2 (from Fig. |4j;), we calculate, (using eqn. 18 of 



iNusser et ail d2006l) and assuming that F — 4/3), that if the bub- 
ble is overpressured by a factor of 2, then the energy imparted 
to the IGM is E IG m 2 = (0.5 - 1.7) x 10 58 erg. If however, 
the bubble is overpressured by a factor of 10, then Eigm 10 = 
(1 - 4) x 10 58 erg. 

To compare this with the radiative losses of the IGM, i.e. the 
X-ray unabsorbed luminosity within the region that the bubbles 
could have affected, a characteristic time-scale needs to be used. 
This time-scale over which this is done is also important, since if 
the energy is injected too slowly, it will be insufficient to counteract 
cooling, and if the energy is released too quickly, it would cause a 
spike in the temp erature of the IGM . 

Following lBirzan et al.l {2004), we define three different time- 
scales, and thus potential ages for the bubble. The first (t c ) is based 
on the time it would take for the bubble to rise to its current loca- 
tion from the BGG at the group's sound speed, if the bubble had 
been inflated by a power ful jet close to the radio core (see for ex- 
ample the simulations of lOmma et al.| [2004). The second estimate 
(tb) is based on the assumption that the bubble is again inflated 
close to the radio core but rises buoyantly rather than being pro- 
pelled due to its production mechanism. The third method is to 
calculate h ow long it would take to re fill the volume displ aced by 
the bubble dMcNamara et alj|2000l a ndlNulsen et alj|2002l) . which 
is denoted t r . We use eqns. 2-4 from lBirzan et al.l d2004h to calcu- 
late time-scales based on the sound-speed (t c ), the buoyant rising 
time of the bubble (tt) and the refilling time for the bubble (t r ) 
respectively. The time-scale and powers calculated are shown in 
Tableland the unabsorbed luminosity of the group within 29 kpc 
is 0.36±0.02 x 10 42 erg s _1 . Thus, it is clear that over the lifetime 
of the bubble, and for the foreseeable future, energy input from the 
bubble will dominate over radiative losses (by 1-2 orders of magni- 
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Figure 5. Radio and optical emission around NGC 741/2. The greyscale shows the //STF814W image described in the text, with grey levels chosen to 
emphasise the common stellar envelope of the system. Contours are from our 1.5-GHz radio image with resolution 3.1 X 1.7 arcsec (major and minor FWHM 
of elliptical restoring Gaussian) at 0.5 X (1, 2, 3 . . . 10, 15, 20 . . .) mjy beam -1 . The three optical galaxies and the bent jets of NGC 742 are labeled. 



Table 3. Energy injection rates based on different methods of bubble motion and bubble formation. 



X 


EiGM 


t c 


h 


tr 






Pt T 


Lx 




10 58 erg 


10 7 yr 


10 7 yr 


10 7 yr 


10 42 ergs" 1 


10 42 ergs" 1 


10 42 ergs" 1 


10 42 erg s" 1 


2 


0.5-1.7 


4.5 


8.6 


6.7 


4-10 


2-6 


4-20 


0.36 ±0.02 


10 


1-4 


4.5 


8.6 


6.7 


7-30 


4-15 


9-30 


0.36±0.02 



t c , P c are the time-scales and powers derived assuming that the bubble was inflated by a powerful jet and rose to its current position at the 
sound speed of the group. 

tb, P;, are the time-scales and powers derived assuming that the bubble was inflated by a more gentle jet and rose to its current position 
buoyantly. 

t r , P r are the time-scales and powers derived from refilling arguments. 



tude), providing more than sufficient energy to counteract cooling 
losses, even if heating is not 100% efficient. 



5 PHYSICAL CONDITIONS IN THE BUBBLE 

If the observed cavities in the X-ray emission correspond to bub- 
bles injected by an AGN outburst, then we can use the X-ray and 
radio data available to constrain the nature of the plasma filling the 
bubbles. Our starting assumption is that the pressure in the bub- 
ble must be at least equal to the external thermal pressure, while 
it cannot be very much greater than that pressure to avoid driving 
a shock into the external medium. In addition, the bubble must be 
light compared to the external medium: it is hard to see how an 
over-dense 'bubble' could have reached its present distance from 
the nucleus of NGC 741. 



5.1 Content of the bubbles 

There are then two basic possibilities for the nature of the bubble 
material. The first is that it is very hot (i. e . Tb ^ > T ext), low-density 
mater ial (for example Mathe ws et al.l 1200 3k nd iMcCarthv et al.l 
2007). The second, and in our view more a priori likely, possibility 
for the nature of the bubble material is that it consists wholly or 
partly of a relativistic plasma, containing at least relativistic elec- 
trons and magnetic field, as found in the lobes of radio galaxies. 

In the first case, the material would simply be a hot thermal 
plasma, and thus, a signature of this gas may be expected in the X- 
ray spectrum of t he bubble. To investigate th is, we used a method 
similar to that of Sanders & Fabian (2006); we extracted a spec- 
trum for the region defined as the cavity in Section [2~71 We fit the 
background subtracted (using a blank-sky background), grouped 
spectrum with an absorbed MeKaL model with a temperature 
fixed to that of the gas immediately around the cavity, with the nor- 
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malisation free to vary. This model provides a reasonable fit to the 
cavity spectrum, but cannot conclusively rule out the presence of 
very hot gas in the cavity. Thus, this component is frozen in the 
fit and a second MeKaL component is added to model the possi- 
ble hot phase with fixed abundance (set to O.3Z0) and temperature, 
leaving only the normalisation of this component free to vary. We 
vary kT by hand between the temperature ranges of 5-40 keV, and 
for each value of kT we obtain a 3 — a upper limit on the normal- 
isation and convert this into the maximum pressure of the hot gas 
at a given temperature. If this pressure is significantly greater than 
the pressure of the external gas, it implies that the bubble would 
not be stable under this configuration (as there may also be an ad- 
ditional pressure component from non-thermal sources). However, 
if the pressure in the bubble from the putative hot component is 
less than the external pressure, this provides a lower limit on the 
temperature of any hot gas component. By varying the temperature 
of the second component until the maximum pressure in the hot 
component exceeds the external pressure, we can set a limit on the 
temperature of any hot component present in the gas. It should be 
emphasised that this is not a 'fit' to the spectrum in the usual sense 
as there are insufficient counts to fit a two-component model. The 
method described merely provides a way of testing for the possible 
presence of a second component. 

As the bubble is in direct mechanical contact with the IGM at a 
range of different radii (due to the fact that it is three dimensional), 
using a single temperature from the deprojected temperature profile 
would not be an accurate reflection of the temperature around the 
cavity. However, using a range of temperatures would not provide 
a satisfactory constraint on the pressure. We thus obtain the pro- 
jected temperature in an annulus at the radius of the centre of the 
cavity, to provide a more representative temperature and pressure 
constraint. Fitting the third annulus from the centre of the Chan- 
dra data with a simple absorbed MeKaL model in the 0.5-5.0 keV 
energy range and fixing the abundance to 0.3Zq, we obtain an ex- 
ternal temperature of 1.41 ±0.1 keV. Assuming that the lower limit 
of the pressure of the gas is approximately 5 x 10~ 12 dyne cm -2 , 
an upper limit for Th can be found. We obtain a limit of 10-20 keV 
for the temperature of any putative hot component; above this tem- 
perature range, the gas is over-pressured with respect to the external 
pressure. Thus, we cannot rule out the presence of hot tenuous gas 
in the bubble. 

In the second case, as discussed in Section [T] the combina- 
tion of electron energy spectrum and magnetic field must be such 
that the bubble produces no detectable radio or inverse-Compton X- 
ray emission. For simplicity we consider a uniform electron spec- 
trum and magnetic field strength throughout the bubble, and take 
the magnetic field to be 'tangle d' so that it has an isotropic mag- 
netic stress tensor dLeahvlll99lh . Let the electron energy spectrum 
be N(E) such that N(E)dE is the number density of electrons 
with energies between E and E + dE, and let the energy density 
in non-radiating particles be some multiple « of the energy density 
in electrons, as is conventional: then the requirement of pressure 
balance gives us 

»E m „ 



Ml 

2^0 



+ (!• 



EN(E)dE = 3p c 



(3) 



where B is the magnetic field strength and E mlrL and i5 max are 
the limits on the electron energy densities to be considered. If we 
further assume some physically reasonable form for the electron 
energy density, such as a power law in energy (N(E) = NoE~ F ), 
then, for a given B, k, E m i n and i? ma x we can solve equation[3]for 
the normalization No of the electron spectrum, and so, effectively, 



for the total number or energy density of the electron (particle) pop- 
ulation. 

Pressure balance has been investigated using this basic for- 
malism for a variety of differe nt types of filled cavitie s, with differ- 
ing results. In powerful FRII dFanaroff & Rilevll 19741) radio galax- 
ies, it has been possible to measure the magnetic field strengths 
in the lobes of some systems by detection of inverse-Compton 
scattering of the cosmic microwave background (CMB): the field 
strengths turn out to be comparable to, though slightly below, the 
equipartition values for an electron/positron (ft = 0) plasma only 
dCroston et al]2005f) . so that B « B eq , where 



2fx 



(l + «) 



EN(E)dE 



(4) 



and the corresponding plasma+magnetic field pressures are, in 
several cases, comparable to the measur ed external pressures 
dHardcastle e t al. 2002s; ICroston et afll2004) . suggesting that ft is 
indeed « in these systems. Even where no external pressure can 
be measured the fact that the measured magnetic field is close to 
the equipartition value for k — has to be a systematic coinci- 
dence if (1 + ft) S> 1 in these systems. On the other hand it has 
been known for many years ( e.g., Morganti et al] ll988l;lFerettietal.l 
ll992l : lHardcastle et alj|l998l : IWorrall & Birkinsha w 2000) that the 
minimum internal pressure (ft = 0, B « B eq ) in the lobes of low- 
power FRI radio galaxies of ten falls orders of m agnitude below the 
measured external pressure: Dunn et al. (2005) have recently com- 
mented on this effect in a sample of cluster-centre objects. In these 
systems it seems that either (1 + ft) ^> 1 (as a change in the fill- 
ing factor of the plasma has a similar effect to an increase in ft, 
we do not consider this separately) or there is a strong departure 
from equipartition. In some cases the lack of detectable inverse- 
Compton emission has been used to argue that it is not possible 
that B <C Bcq- In the following section we discuss the constraints 
that the synchrotron and inverse-Compton limits on the NGC 741 
cavity place on the nature of the cavity-filling plasma. 



5.2 Constraints on the nature of the plasma 

The key free parameters of eq. [3] are B, ft, i? m i n and £ max , and 
p. We cannot constrain all of these from our observations. How- 
ever, it seems unlikely that the plasma is magnetically dominated, 
and so one constraint that we can impose is that the magnetic field 
should be less than or equal to the equipartition value, B < B cq . 
The limits on observed synchrotron and inverse-Compton emissiv- 
ity place additional constraints on the parameters of the plasma: 
synchrotron emissivity at a given frequency constrains a com- 
bination of No, B and E m&K , while inverse-Compton emissiv- 
ity constrains No and to some extent -E max only. It is easiest to 
calculat e these constrain t s num erically, and we do this using the 
code of Hardcas tle et al. I fl998h . assuming that the dominant pho- 
ton field for inverse-Compton scattering is the cosmic microwave 
background radiation (i.e. ignoring the negligible synchrotron self- 
Compton component). From the radio data in hand we set an up- 
per limit on the bubble flux density at 1.4 GHz of 7.4 mjy and at 
330 MHz (Jetha et al. in prep) of 20 mJy. These limits were ob- 
tained from VLA observations of the source, with the 1 .4 GHz data 
being the same as used in Section [3] and the 330 MHz from new 
observations (Jetha et al, in prep.). For each observation, since part 
of the radio source NGC 742 obscures the cavity, the non-detection 
limits of the cavity were taken to be three times the off-source rms 
noise multiplied by the square root of the number of beams in the 
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cavity at each frequency. From the X-ray observations, we impose 
an upper limit on the inverse-Compton flux density of 5.8 nJy at 1 
keV. 

We fix the power-law index p at a relatively flat va lue (we 
choose p = 2.1 based on the results o f lYoung et alj|2005l for FRI 
radio sources) which reflects expectations from particle accelera- 
tion models: although there should almost certainly be a tail of par- 
ticles with a steeper value of p due to spectral ageing (as reflected 
in the steep spectral index of some observed radio lobe systems) it 
is very unlikely that it is correct to extrapolate this spectral index 
back to low energies, where the electron loss time-scales are very 
long, and the steep-spectrum tail makes a negligible contribution to 
the overall energy budget of the electrons. We conservatively take 
the minimum electron energy to be E m i n — m c c 2 , i.e. 7 m i n = 1, 
where 7 is the Lorentz factor of the electron). For a given value of 
k, we can investigate the constraints on i5 max (7max^ c c 2 ) and B 
given by the observations on the assumption of pressure balance as 
in equation [3] 

The results are plotted in Fig. |6]for two representative values 
of k as defined above. We chose two values, k = (i.e. relativistic 
electrons/positrons and field give the only contribution to the inter- 
nal pressure of the plasma, as in FRII sources) and k = 10 (i.e. 
there is a population of non-radiating particles, such as protons, 
with energy density exceeding that of the electrons by an order of 
magnitude). We carry out the calculations using the best estimate of 
the external pressure (7.7 x 10~ 12 dyne cm -2 ) and the range im- 
posed by the uncertainties in temperature, density and cavity posi- 
tion (3-10.4X 10~ 12 dyne cm -2 ). It can be seen that the radio and 
X-ray data already impose quite severe constraints on the possible 
electron energy spectrum if the pressure is close to its best-fitting 
value. For n = the combination of inverse-Compton and syn- 
chrotron constraints requires 7 max < 1000 for any plausible value 
of B. Even for n = 10 the magnetic field strengths must be very 
low compared to their equipartition values for 7 max to be greater 
than a few x 10 3 : this compares to 7 max ^2 x 10 4 for the lobes of 
typical radio sources with detections at GHz frequencies (assuming 
fields close to equipartition). For the values of n investigated (0 and 
10), the bubble can be filled with a plasma with 7, nax -* 2 x 10 3 only 
for the lowest external pressures and the lowest internal magnetic 
fields. It is clearly possible to choose a value of k such that we have 
no constraint on the properties of the electron population - this cor- 
responds to a situation in which the energetics of the bubble plasma 
are entirely dominated by an 'invisible' component. However, it is 
interesting to consider the constraints on the electron population for 
a situation in which (1 + k) ~ 1. For these low values of n, how 
can we obtain the required low values of 7 max ? 

Radiative losses alone (to synchrotron and inverse-Compton 
emission) would require a long time to reduce an initially high 
7 max to values ~ 1000. The maximum energy in an electron popu- 
lat ion that initially ex tended to high energies after a time t is given 
bv lPacholczvkl fl97(t > 



where m e is the mass of the electron, gt is the Thomson cross- 
section and U is the total energy density in magnetic field and mi- 
crowave background photons (i.e. the total energy density impli- 
cated in loss processes). For equipartition magnetic field strengths 
synchrotron losses dominate, U = 3p cx t/2 and the age of the 
bubble is > 2 x 10 8 years (though note that this calculation 
takes no account of the fact that the magnetic field may have 
been higher at earlier times). For magnetic field strengths much 



less than equipartition inverse-Compton losses dominate, U = 
87r 5 (fcTcMB) 4 /15(/ic) 3 , and the limit on the age of the bubble is 
closer to > 2 x 10 9 years, where t is the time since any high- 
energy particle acceleration took place in the source. These time- 
scales are considerably longer than the time-scales for the bubble 
to rise buoyantly to its current position (8.6 x 10 7 years, Section[4]l. 
7 max may also be reduced by adiabatic expansion, which gives 
7 max oc y~( r_1 ) 5 where T is the adiabatic index, 4/3. To take 
a simple example, we could suppose that the original bubble had 
7 max = 2 x 10 4 as quoted above and has expanded adiabatically. 
Then the volume must have changed by a factor 8000 during the 
expansion; i.e. the initial radius of the bubble was around 1 kpc. 
But this would require the initial bubble to have been over pres- 
sured by a factor 1.6 x 10 if it is in pressure balance now, which 
seems unlikely. As shown by Alexander (1987), more modest ex- 
pansion factors coupled with radiative losses do not greatly change 
the inferred ages. 

We can therefore conclude that it seems hard for the bubble 
to have evolved from a 'normal' region of radio-emitting plasma - 
like the plumes of NGC 742 whose radio emission we see in the 
system at present - if the lobe external pressure is close to the best- 
fitting value. The present X-ray data do not allow us to constrain 
the pressure around the cavity well enough to rule this possibility 
out completely, though, and a deeper X-ray observation is required 
to allow us to draw firm conclusions: if the pressure is at the low 
end of the range we estimate, and the internal pressure of the cavity 
is dominated by non-radiating relativistic protons with a very low 
magnetic field, then the available constraints are much weaker. In a 
future paper we will apply the technique used above to constrain the 
nature of the cavity-filling plasma in other 'ghost cavity' systems. 



6 CONCLUSIONS 

We have examined the energetics of the galaxy group associated 
with NGC 741. We find that if the cavity seen in the X-ray data was 
originally inflated by a radio-loud AGN associated with the bright- 
est galaxy, NGC 741, then there is sufficient energy to counteract 
cooling. This interpretation is consistent with the fact that NGC 741 
clearly has an active nucleus at the present time, and this could have 
been more active in the past. 

We have investigated the nature of the cavity-filling plasma. 
We cannot rule out the possibility that the cavity is filled with a very 
hot thermal plasma (kT > 10 keV) or with a relativistic plasma 
that is energetically completely dominated by non-radiating parti- 
cles. However, it is not at all clear that current models for AGN 
inflation of cavities can produce cavities of this kind, particularly 
when we take into account the relatively small distance (30 kpc) 
between the cavity and the centre of the host group. If on the other 
hand there is a substantial contribution to the energy density in the 
cavity from a plasma that contains relativistic electrons and mag- 
netic field, as in the active lobes of radio galaxies (such as NGC 742 
in the same group) then we find that the limits placed by the non- 
detection of radio synchrotron and X-ray inverse-Compton radia- 
tion severely constrain the maximum electron Lorentz factor. In this 
scenario, we find that it is very hard for the cavity to have evolved 
to its current state from a 'dead' radio lobe by a combination of the 
standard processes expected to be operating (synchrotron, inverse- 
Compton, and adiabatic expansion losses); for radiative losses only 
the lobe would have to be older ( > 2 x 10 s years) than any plausi- 
ble bubble inflation time (between 5 x 10 7 and 9 x 10 7 years). The 
extreme parameters required for the NGC 741 cavity, and for other 
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Figure 6. Constraints imposed by the radio and X-ray limits on the physical conditions of any relativistic plasma in the bubble. The white regions indicate 
those regions of the parameter space allowed given the constraints, and the greyscale indicates the number of constraints that are violated, ranging from 1 
(light grey) to 4 (both radio limits, the IC limit, and the constraint that the plasma should not be magnetically dominated) (dark grey). From top to bottom, 
the panels show the constraints for external pressures of 10.4 X 10~ 12 , 7.7 X 10~ 12 , and 3.0 X 10 -12 dyne cm~ 2 . The left-hand panels show the k = 
constraints while the right-hand ones show k = 10. Note the different magnetic field and 7max scales of the three sets of plots. 
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'ghost cavities' whose physical parameters have been investigated 
(e.g. Dunn et al 2005) and the fact that cavities are found at a wide 
range of radii (Birzan et al 2004) may indicate that inflation by the 
lobes of a standard radio-loud AGN is not a viable explanation for 
all ghost cavities. In the case of NGC 741, we are left with one 
clearly outstanding question: if the cavity in NGC 741 is due to an 
AGN-inflated bubble, in what form was the energy injected into the 
IGM, and in what form is it now? 
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